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Abstract 

The pp interactions taking place in the cosmos around us are a source of the astrophysical 
neutrinos of all the three flavors. In these interactions, the electron and the muon neutrinos 
mainly come from the production and the decay of the tt^ mesons, whereas the tau neutrinos 
mainly come from the production and the decay of the mesons. We estimate the three intrinsic 
neutrino flavor ratios for 1 GeV < E < 10 12 GeV in the pp interactions and found them to be 1 
: 2 : 3 x 10~ 5 . We study the effects of neutrino oscillations on these intrinsic ratios. We point 
out that the three ratios become 1 : 1 : 1 if L(pc) / E(GeV) > 10~ 10 in the presence of neutrino 
oscillations, where L is the distance to the astrophysical neutrino source in units of parsecs. 
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I. INTRODUCTION 



Neutrino astronomy holds a great promise to explore the interiors of the dense astrophys- 
ical systems, which is not possible by any other existing means such as through the study 
of the cosmic-rays and/or the gamma-rays . It is mainly because once the neutrinos 

are produced in the distant cosmos, they are essentially unobstructed by the intervening 
background matter mainly owing to their weak interaction cross section. On the other hand, 
the cosmic-rays (being the charged particles) and the gamma-rays are either deflected or 
even considerably absorbed by the same intervening background matter, depending upon 
the energy j^. 

The presence of the proton component in the observed cosmic-ray flux for the entire 
energy range (1 GeV < E p < 10 12 GeV), may already be signaling the anticipated existence 
of the neutrino astronomy. The search of neutrinos from the cosmos may help to find a 
unified explanation for the common origin of the ultra-high energy cosmic-rays and the 
high energy gamma-rays. The absolute levels of the astrophysical neutrino productions 
are determined/affected by the intervening background matter P, |(|, the relative levels 
nevertheless do not. These thus remain important observables for the forthcoming detailed 
astrophysical neutrino searches. It is thus important to investigate the relative neutrino 
production levels in an astrophysical neutrino source; as well as the changes that may occur 
in the relative composition of the neutrino flavors during their propagation to us. 

Given the recent empirical evidences of the neutrino oscillations p|, it is timely to 
perform a reference study for the effects of the neutrino oscillations on the mixed intrinsic 
astrophysical neutrino flavor ratios during their propagation. Further motivation is provided 
by the recent developments both in providing examples to use the astrophysical neutrinos 
not only to study the cosmos around us ^, a, 10, [ll|, but also to explore the properties 
of the neutrinos itself (including those suggested beyond the Standard Model of Particle 
Physics) 0. 

We consider the pp interactions as a source of the intrinsic astrophysical neutrino pro- 
duction. The first p represent the cosmic-ray flux produced inside the source, whereas the 
second p represent the medium contribution in the source. Commonly cited examples of the 
astrophysical systems where the pp interactions play a role include the nearby astrophysical 
sources such as the earth atmosphere and the galactic plane/center region as well as other 
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sources inside our galaxy. The more distant suggested astrophysical neutrino sources include 
the Active Galactic Nuclei (AGNs) 3, and the sites of the Gamma-Ray Bursts (GRBs) 

3- 

We estimate the three intrinsic neutrino flavor ratios as a function of the neutrino energy 
and study the neutrino oscillation effects on these for neutrino energy ranging between 1 GeV 
and 10 12 GeV. The energy dependence of the three intrinsic ratios has not been studied 
previously. We mainly investigate the particle physics aspects of the three ratios. A purpose 
of the present study is to provide a firm basis for the relevance of the neutrino oscillation 
effects for the forthcoming searches of the three astrophysical neutrino flavor ratios. As in 
order to determine the astrophysical source characteristics as precisely and as completely as 
possible via neutrinos, one needs to know the relevance of the distance to the source L in 
the presence of the neutrino oscillations for the given neutrino energy E, since the neutrino 
oscillations redistribute the intrinsic neutrino flavor ratios depending upon the value of the 
ratio L/E. 

This paper is organized as follows. In Section II, we estimate the three intrinsic astrophys- 
ical neutrino flavor ratios in some detail mainly within the framework of the Quark-Gluon 
String Model (QGSM). In Section III, we study the effects of the neutrino oscillations on 
these. We identify the range of the L/E values where the commonly considered assumption 
of the averaging of neutrino oscillation probabilities may hold. In Section IV, we briefly 
summarize the present status of the relevant detection strategies as well as the detector 
developments to search for the astrophysical neutrino flavor ratios. In Section V, we present 
our conclusions. 



II. THE INTRINSIC ASTROPHYSICAL NEUTRINO FLAVOR RATIOS 

Let us define the three intrinsic neutrinos flavor ratios as the R®/ e , the R^u, arid the 
R° r/e . Here R° r/e = F° r (E)/F^(E) with F? r (E) = dN°JdE : for instance. Clearly, R° e/e = 1. 
This ratio provides normalization for the other two ratios. The previous estimates for the 
R° T j e ratio are between > 10~ 4 ^| and > 10 -5 Q]. However, the energy dependence of the 
three intrinsic neutrino flavor ratios was not studied. 
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We use the following formula for computing the astrophysical neutrino flux spectrum 

= ^ = dn p ^ dE p ME P ) , (1) 



where E is the neutrino energy. The cosmic-ray flux spectrum, <ft p (E p ), is given by 17 1 

<j) p (E p ) = A(E p /GeV) s cHrW^GeV" 1 , (2) 

where A = 0.8 and 5 = —2.75. We use the above cosmic-ray flux spectrum for 1 GeV 
< E v < 10 12 GeV. We are aware that the above cosmic-ray flux spectrum differs by few 

~ n 

percent for E p < 10 GeV as compared to the more recent compilation [18||. However, this is 
of not much concern for our present study since we are primarily interested in studying the 
astrophysical neutrino flavor ratios here. Moreover, later, we also study the effects 
of varying the exponent 5. In the above simplified picture, it is assumed that all the hadrons 
and the relevant leptons decay before interacting with the medium of the astrophysical 
neutrino source. 

The representative proton number density inside the source is taken to be n p = 1 cm -3 , 
and the representative distance d inside the source is taken to be ~ 10 kpc, where 1 pc ~ 
3 x 10 18 cm. These values are to merely represent the reference absolute levels for the three 
neutrino fluxes. As stated earlier, our main concern in this paper is to study their ratios 
and their energy dependence defined earlier which are obviously independent of the product 
dn p . It is clear that the task of computing the dN^/dE in Eq. (0) essentially relies on the 
evaluation of the differential cross section da/dE in the pp interactions. 

In this work, we shall consider only the most dominant production channels for each 
neutrino flavor generation as the representative examples. Namely, the 7T ± meson for the 
electron and the muon neutrinos and the Df meson for the tau neutrinos. We employ the 
Quark-Gluon String Model (QGSM) to calculate the production distributions of the above 
mesons. The QGSM approach is non-perturbative and is based on the st ring fragmentation. 



It contains a number of parameters determined by the experiments |ial2fl|2ll|22j,:23]. 
The production cross section of the hadron h in the QGSM is given by 
da h (s,x) 1 



^(s) ( j ) h n (s,x)+a DD (s) ( j ) h (s,x) 



n=l 



(3) 



where x = 2p|/-y/s and x± = 2^J(m 2 h + pf)/s. The p| (p^) is the parallel (perpendicular) 
momentum of the secondary hadron h in the center of mass frame. All the related formulas 
are provided in the Appendix. 
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We have used the VEGAS multi-dimensional Monte Carlo integration program 
to obtain the energy distribution of the neutrino flux by selecting the events falling inside 
the considered energy segment through a series of boosts to the laboratory frame. We have 
simplified the secondary decays of the charged leptons, the \x and the r, as effectively the 
two body ones, that is, a neutrino plus a particle with varying invariant mass. For instance, 
we assume that the r lepton decays into a v T and a particle Y with the mass my satisfying 
0.1 GeV < my < m T — 0.1 GeV. In our setting s ~ 2m p E p . We take the Dg — > t + v t 
branching ratio as ~ 0.064 |26j ]. 

We have calculated the following processes in the QGSM: 

pp -> D±(-> tv t )+X 

> v T + y. (4) 

The tau neutrino production via the Z)^ meson can also be calculated in the perturbative 
Quantum Chromo Dynamics (pQCD) through pp — > cc — > -D^ + X — > z/ T + K . However, the 
electron and the muon neutrino production through the 7r ± meson can not be handled in 
the pQCD. The reason is as follows. There are many quark-level sub-processes (the the 
u- as well as the s-channel) for the tt^ production in the pp interactions, differently from 
the Df production which has only the cc pair production channel. Without some cuts such 
as on the or some cuts on the factorization scale etc., the light quark productions may 
blow up as they approach the non-perturbative region. 

In the pQCD calculations, we use the leading-order results of the parton sub processes 
QQj 99 ~^ QQ where Q = c for pp — > (cc — > Df — >) v T Y [and Q = t for pp — >• ti — >• vX\. 
We use a K factor, K — 2, to account for the NLO corrections. For the parton distribution 



functions, we use the CTEQv6 |22ll2S|. We use m c = 1.35 GeV, m t = 175 GeV, a s with 
a s (M|) = 0.118 and Q 2 = s/4 as a factorization scale. For the v T production through the 
cc, we use the Peterson fragmentation function with e ~ 0.029 for fragmentation of c or c 
into the Dr- meson 

HQ- 
Fig. H shows the three intrinsic neutrino flavor fluxes as a function of the neutrino energy 

E. We plot diV°/d(log 10 £) in units of cm 2 s x sr 1 as a function of the E. Since the tau 

neutrino production via the D§ can be dealt with in both the pQCD and the QGSM, the 



pQCD result for the tau neutrinos is also presented for comparison. We note that there is 
a relatively large discrepancy at higher energy (E > 10 9 GeV) between the two approaches 
for the v T production. 

To consider an example of the process that may produce the three neutrino flavors without 
hadronizing in the pp interactions, we have studied the pp — > tt channel. Here, the neutrino 
energy distributions are not effected by the hadronization process. The neutrino production 
in the pp — > tt — > uX is reliably calculable in the pQCD. Here, we consider the direct decays 
of t into each lepton (with ~ 10% branching ratio) and include the secondary decays of 
massive leptons (/i and r). Note that all the lepton masses are very small as compared to 
the top quark mass, as a result all the three neutrino distributions are of the same orders 
of magnitude. These results are also shown in the Fig. It is clear from the figure that 
the three neutrino fluxes are comparable for this process. Compared to the 7r ± and the 
.Dg results, the tt contributions are negligible, less than a factor of at least ~ 10~ 3 over the 
whole range of the considered energies. 

It can also be seen from the figure that the production of the tt^, the Df, and their 
branching ratios of relevant leptonic decays are more important than the effects of the 
lepton masses. Our results are in good agreement for the v e and the production with 
those given in the Ref. [3jj using the PYTHIA, whereas for the v T production, our results 
are in good agreement with those given in the Ref. 

Fig. 121 shows the three intrinsic neutrino flavor ratios defined earlier, as a function of the 
neutrino energy E. We note that an energy independent relative flux hierarchy among the 
three intrinsic neutrino flavor ratios persisted even at the highest considered energy, namely 
R e/e : R l/e '■ K/e oc 1 : 2 : 3 x 1(T 5 for 1 GeV < E < 10 12 GeV. Our results are however 
subject to the uncertainties of extrapolating the parameters of the hadron/quark production 
models especially for E > 10 6 GeV. This corresponds to center-of-mass energy of y/s > 10 3 
GeV. To our knowledge, the QGSM parameters are fitted up to y/s ~ 540 GeV using the 
SPS collider data for the light mesons, whereas for the charmed mesons the comparison is 
available using the y/s ~ 630 GeV data. 

In Fig. El we show the energy dependence of the two intrinsic ratios of the astrophysical 
neutrino fluxes, the R^/ e and the R®/ e , by varying the exponent 5 of the cosmic-ray flux 
spectrum, <fi p (E p ) [see Eq. (J2J)], in the range —1.75 < 5 < —3.75 . It is so because the cosmic- 
ray flux spectrum in an astrophysical source is expected to be harder than the observed one 
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locally 17(. We note that the ratio R° T / e changes from ~ 3 x to ~ 8 x 10~ 5 , when the 
5 changes from —2.75 to —1.75 for 1 GeV < E < 10 12 GeV. However, except the above 
(slight) change, in general, we note from the figure that the two intrinsic neutrino flavor 
ratios are essentially stable w.r.t. the 5 variation. 

For an example of possible energy dependence in the the three astrophysical neutrino 
flavor ratios coming from astrophysical reasonings, see 3^ |. 



III. THE OSCILLATED ASTROPHYSICAL NEUTRINO FLAVOR RATIOS 

In this Section, we shall perform a three neutrino oscillation analysis of the three intrinsic 
neutrino flavor ratios estimated in the previous Section. In the context of the three neutrino 
flavors, there are 6 independent neutrino mixing parameters. The matter effects are found 
to be negligible for the entire range of the neutrino mixing parameters and the E values 
under discussion 0, 3fi| . 

In this analysis, we do not make the assumption of averaging over the neutrino oscillation 
probabilities. The neutrino oscillation effects for the astrophysical neutrinos, using the 
averaged oscillation probability expressions, were studied in some detail in Ref. Here, 
instead, we shall determine the L/E range that may be relevant for the averaging. 

We start with the connection U between the flavor | u a ) and the mass | z/j) eigenstates of 
the neutrinos, namely 

3 

I "a) = ^ Uai 1 Vi), (5) 
1=1 

where a = e, fi, r. In the context of the three neutrinos, U is called the Maki Nakagawa 
Sakita (MNS) mixing matrix [3j|. Under the assumption that the CP violating phase 
5cp = 0, the 3x3 MNS mixing matrix U in the standard parameterization connecting the 
neutrino mass and the flavor eigenstates reads 



/ C12C13 S12C13 S13 1 



U 



(6) 



— S 12C23 ~ C12S23S13 C12C23 — Si2S23 s 13 S23C13 
y S12S23 ~ c 12 c 23 s 13 — C12S23 — S12C23S13 C23C13 J 

where Cy = cos 6^- and sy = sin 6^-. The presently available empirical constraints for the 
various neutrino mixing parameters, in the context of the three neutrino mixing, give the 
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following best fit values: 12 = 33.2°, 9 23 = 45.0°, 9 13 = 0.0°, 5m{ 2 = 7.9 x 1CT 5 eV 2 and 



0. 

Using Eq. ©, the neutrino oscillation probability formula is |3£ 



5m 2 23 = 2.1 x lCr 3 el/ 2 



P(z/ a - ^; L, E) = P aP {L, E) = J2 + ^ U aj U pi U aj U Pj cos , (7) 

i=l i# \ L ijJ 

where (3 = e,fi,r and L^- ~ AE/Smfj is the neutrino oscillation length. The L in Eq. (|7J) is 
the neutrino flight length. 

The neutrino flux ratios R a /p(L, E), arriving at the detector, in the presence of neutrino 
oscillations are estimated using the relation 

F Va {L,E) = Y, P ^ L ' E ) F ^ ( 8 ) 

P 



and the definition R a /p = F Ua (E) / F U0 (E) , so that 

K a /f3{L,E)- - — T ^ | - lT fNd0 1 j~) (j tp\ do ? (9 J 



PgejL, E) + P Q „(L, E)i?; /e + P ar (L, P)P° /e 

P^L, P) + P^(L, P)Pj /e + P^L, P)P° /e ' 



where the P^ are taken according to the discussion in the previous Section. The P a p(L, E) 
is obtainable using Eq. (JJJ). The unitarity conditions such as 1 — P ee (L,E) = P efl (L,E) + 
P eT (L, E) are implemented at each L and E at which these are evaluated. 

Fig. |3] shows the three oscillated neutrino flavor ratios as a function of the ratio L/E. 
Here, we have used the fact that the three intrinsic ratios are essentially independent of the 
L as well as the E. Note that the ratio Rr/^ behaves like the ratio R T / e . The ratio R e / e is 
not plotted as it is unaffected by the neutrinos oscillations. From the figure, it is clear that 
the averaging may be a good approximation if L(pc)/P(GeV) > 10 -10 , namely 

1 : 2 : 3 x 1(T 5 — ► v osc and if L(pc)/E(GeV) > 1CT 10 — ►1:1:1. (10) 

Fig. which is our main result, shows the region in the L versus E plane where the 
averaging of the neutrino oscillations probabilities may be assumed, under the above crite- 
rion. For instance, if E ~ 10 4 GeV, and if the distance to the source is > 10 12 cm then the 
incoming astrophysical neutrino flux should be essentially an equal admixture of the three 
neutrino flavors. 
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IV. PROSPECTS FOR ASTROPHYSICAL NEUTRINO FLAVOR IDENTIFICA- 
TION 



In this Section, we shall briefly summarize the presently envisaged detection strategies 
and the detector configurations for the possible astrophysical neutrino flavor identification 
for 1 GeV < E < 10 12 GeV mainly in the context of the Cherenkov radiation detection. For 
a discussion of other alternative detection strategies such as using the radio and the acoustic 
signals, see Ref. HQ. 

From the reference estimates presented in Section II, it follows that the flux of the intrinsic 
astrophysical tau neutrino flavor is considerably suppressed for 1 GeV < E < 10 12 GeV rela- 
tive to that of the electron and muon neutrino flavor. According to the discussion in Section 
III, the neutrino oscillations populate the astrophysical tau neutrino flux comparable to the 
astrophysical electron and the muon neutrino fluxes, provided the ratio L/E is in a certain 
range. A signature of the neutrino oscillations in the mixed astrophysical neutrino flux thus 
shall be the identification of the astrophysical tau neutrino flavor among the 
other two. A commonly used essential ingredient in this context is to make a comparative 
use of the characteristic energy dependent astrophysical tau neutrino induced tau lepton 
decay and/or interaction length scale relative to the relevant electron and the muon length 



scales 
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A. Current detection strategies 

The astrophysical neutrino detection can be achieved in the neutrino nucleon/electron 
interactions QQ. These interactions may occur near or inside the detector. The charged 
leptons, the (air) showers, and the associated radiations such as the radio and/or acoustic 
signals are the measurable quantities. The detectors are optimized for their performance in 
discriminating the three neutrino flavors in certain energy intervals 46]. In addition to re- 
constructing the astrophysical neutrino flavor ratios from an individual detector, comparing 
the data from the various detectors shall also lead to the possibility of identification of the 
three neutrino flavor ratios. For a recent discussion of the event rates in some representative 
astrophysical neutrino flux models, see Ref. j^. The astrophysical neutrinos arrive at an 
earth based detector in the three general directions. 
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The downward going astrophysical neutrinos do not traverse any significant cord length 
inside the earth to reach the detector. In fact, for large earth surface shower detectors, the 
neutrino nucleon interactions take place in the earth atmosphere and the resulting shower 
(or part of it) is observed by the detectors. Several studies were performed to identify the 
detector specifications and/or the energy ranges in which a specific detector configuration 
is/can be optimized to disentangle the astrophysical tau neutrino flavor from the astrophys- 
ical electron and /or muon neutrino flavor through the double shower or single shower event 



topologies 



HQ 
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The upward going neutrinos traverse a large earth cord before they reach the detector. 
At energies E > 5 x 10 4 GeV, the charged current neutrino nucleon interaction length 
is smaller than the earth diameter. As a result significant neutrino flavor dependent 
absorption takes place for E > Eq. For a discussion of the upward going; tau neutrino 
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53]. 



flavor behavior versus the upward going muon neutrino flavor, see Ref. 

It might also be possible to search for the air showers induced by the incoming (quasi- 
horizontal) astrophysical neutrinos, in case the neutrinos happen to have just one interaction 
inside the earth. This strategy is referred to as the earth skimming £4( . In some configura- 
tions, the astrophysical neutrinos may interact just below (within 5° ~ 10° of) the detector 



horizon 
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]. The air showers produced by the earth skimming astro- 
physical neutrinos can also be searched by the forthcoming large scale balloon/space based 
detectors as well 



B. Present and the forthcoming detectors 



Presently operating detectors searching for the astrophysical neutrinos include the 
Antarctic Muon and Neutrino Detector Array (AMANDA), its proposed extension, the 
Ice cube 0,0,0, and the lake Baikal detector These detectors use the ice and the 
water as detection medium respectively, and are sensitive to all the three neutrino flavors 
essentially for 10 3 GeV < E < 10 6 GeV and mainly search for the upward going (and hor- 
izontal) neutrinos. Other under construction large scale detectors include the Astron omy 
with a Neutrino Telescope and Abyss environmental REsearch (ANTARES) project [67]. 
For E < 10 3 GeV, the Superkamiokande and the upcoming one Mega ton class of detectors 
shall be sensitive to the three neutrino flavors 681. 
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Several attempts are underway to implement the earth skimming strategy using the 
mountains as target for the neutrino nucleon/electron interactions, such as the concept study 
carried out by the Neutrino Telescope (NuTel) collaboration 69j. This class of detectors 
are/shall be essentially sensitive for 10 6 GeV < E < 10 8 GeV, mainly for the tau neutrino 
flavor. Earth/air/sea skimming tau neutrino air shower search for more wider energy range, 



namely for 10 GeV < E < 10 GeV, is also recently suggested 
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Zl|. 



The under construction large surface array detectors such as the Pierre Auger (PA) 



observatory [72|, |73|, and the Telescope Array (TA) experiment shall also be sensitive to 
all the three neutrino flavors for 10 7 GeV < E < 10 11 GeV Q Orbiting Wide-field 
Light collector space based mission (OWL) shall be sensitive to the z/ e (z/ e ) flavor for 10 10 
GeV < E < 10 12 GeV as it shall search for the atmospheric fluorescent trail using earth 

n 

atmosphere as the detection medium in the neutrino nucleon interactions |75[ . The Extreme 
Universe Space Observatory (EUSO) shall be sensitive to all the three neutrino flavors by 
detecting the fluorescent and the Cherenkov light produced in the air showers generated by 
the neutrino nucleon interaction occurring in the earth atmosphere for 10 10 GeV < E < 10 12 
GeV Q 

The detectors based on the alternative techniques are also taking data. These include 
the Radio Cherenkov Experiment (RICE) that is sensitive to the v e [pe) flavor based on the 
anticipated radio-wavelength Cherenkov radiation detection that shall be produced by the 
neutrino nucleon interactions in the polar ice for 10 7 GeV < E < 10 12 GeV 
high altitude balloon based Antarctic Impulsive transient Antenna (ANITA) experiment 
shall search for the ice skimming z/ e (z/ e ) flavor induced coherent radio signals for E > 10 9 
GeV For E > 10 11 GeV, upper limits are also provided by the the Goldstone Lunar 
Ultra-high energy neutrino Experiment (GLUE) mainly for the v e [y e ) flavor, based on the 
similar detection technique |80j]. For a summary of upper limits based on the alternative 



detection methods for E > 10 9 GeV, see |81|. 



V. CONCLUSIONS 



We have performed a reference estimate of the three intrinsic astrophysical neutrino flavor 
ratios for the neutrino energy ranging between 1 GeV and 10 12 GeV in the pp interactions 
mainly within the framework of the Quark-Gluon String Model (QGSM). We have taken into 
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account the tt meson production for the generation of the v e and the neutrino flavors, and 
the meson production for the generation of the v T neutrino flavor. We have also studied 
the tt production using the perturbative Quantum Chromo Dynamics (pQCD) as an example 
of the process for the generation of the three neutrino flavors without hadronization in the 
pp interactions. The neutrino generation from the latter channel is found to be suppressed 
for the entire considered energy range. 

We have only taken into account the proton component in the observed cosmic-ray flux 
spectrum. We have also studied the variation of the cosmic-ray flux spectrum exponent 5 
on the three ratios and found that the three intrinsic ratios are essentially independent of 
the exponent for —1.75 < 5 < —3.75. 

The three astrophysical neutrino flavor ratios are essentially independent of energy 
and are 1:2: 3xl0~ 5 for 1 GeV < E < 10 12 GeV. Namely, the relative intrinsic neutrino 
flux hierarchy stays the same for the entire considered energy range. Therefore, the intrinsic 
astrophysical tau neutrino flavor is relatively suppressed for the entire considered energy 
range. Our considered energy range covers the entire E range of the observed cosmic-ray 
flux. 

We have studied the effects of the neutrino oscillations in the three neutrino 
flavor framework on the three intrinsic ratios, using the recent best fit values of the neu- 
trino mixing parameters. The neutrino oscillation effects depend upon the distance to the 
astrophysical source L for the given neutrino energy E. For L(pc)/£'(GeV) > 1CT 10 , the 
averaging of the neutrino oscillation probabilities may be assumed, where the L is in the 
units of parsecs. Our present estimate is intended to provide a firm basis for the relevance 
of the neutrino oscillations effects for the forthcoming search of the astrophysical 
neutrino flavor ratios by the various detectors. 
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APPENDIX A: FORMULAS FOR THE QGSM 



In Eq. (J3J), the functions a pp and <j>^(s,x) are defined as follows: 



KM 



a h [F q h q (x + ,n)F q h (x^n) + F q h (x + ,n)F q h q (x^n) 
+2(n-l)F^ a (x + ,n)F^ ea (x.,n)] for n > 1, 



where x± = (a/x 2 + ± x)/2 and 



F q h (x,n) = 2 -j\ Xl f^ n )G^(^ +1 f«( Xl ,n) G h d (j^J 



F q h q (x,n) 



2 

3 . .. 

| J d Xl f?( Xl , n) G h ud + dx, f^(x u n) G 



h 

uu 



4 + 2<L 



dzi f^ a {xx,n) 



X 

X 
X\ 



(Al) 
(A2) 



ct \-\ Ct [ -) 

XxJ 



+ / dxt /^ ea (xi,n) 

J X 



X 



x 



+6 S / dxi f; sea (xi,n) 



G h A-\+G h A- 



(A3) 



In the above equations, fp(x, n)'s are the distribution functions describing the n— Pomeron 
distribution functions of quarks or diquarks (i = u, d, uu, ...) with a fraction of energy x from 
the proton, and G\ (z)'s are the fragmentation functions of the quark or diquark chain into 
a hadron h which carries a fraction z of its energy. 
The list of the /*(x,n) is as follows 

Y{l + n-2a N ) 



/; u (x,n) 



f" sea (x,n) 



T(l - cm) r(a fl - 2«at + n) 

r(2 + n - 2a N ) 
T(l — a ft) T(a R — 2a N + n + 1) 



x <£-°ft — x ^ a R~ 2a N+{n-l) 



x x _Qii (1 - X 



,a R -2a N +n 



T(2 + n-2a N ) 



T(-a R + n) T(a R - 2a N + 1) 
r(l + n- 2«7v) 



X X 



OA-20JV+1 ^ _ ^ -aR+(n— 1) 



X X 



r(-a fl + n) T(a R - 2a N + 2) 

T(l + n + 2«r — 2«at — 2a«^) 
T(l — a^) r(2aij; — 2a7v + n — a^) 
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(1-x) 



-a s +(n-l) 



x x - "^ (1 — x) 2 " ii ~~ 2ajv+ ( n ~ 1 -'~ Q ^ 



where T is the usual Gamma function. 
The list of the G? s (z) is given by 

G^{z) = (i- z ) A -«*+ 2 ( 1 -«*), 

= <??'», 
Gf°(z) = (1 - z) x ~ a ^ x (1 + ai z 2 ), 

= Gf{z). 

The list of the Gf{z) is given by 

Gf(z) = {l-z)-^ 

= Gf(z) = Gf(z) = Gl(z), 

Gf(z) = (l-z)- aR+x+1 

= Gt(z) = Gl(z) = Gf(z), 

Gl~( Z ) = (1 - z )°*-2°"+A+l 

In the above, the input parameters are as follows: 

an = 0.5 , cin = —0.5 , 0(0 = , — —2.18 , A = 0.5 , a\ — 5 
5 S = 0.25 (0), if strange sea contribution is turned on (off), 
a D ° = 0.0007 , a w = 0.44 . 

The function <J pp (s) and a DD (s) are given by the following formulas: 

*(0 = £ (i-«p<-*)E£) • 



ODD® = ^° P \f(*m - /(*)] , 



where 
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^ ™ z Jo x 
The best fit parameter values are as follows : 

(i) for y/s < 10 3 GeV 

7 P = 3.64 (GeV)" 2 , R 2 = 3.56 (GeV)~ 2 , a' p = 0.25 (GeV)" 2 , 
C= 1.5 , A = 0.07 . 

(ii) for yfi > 10 3 GeV 

7p = 1.77 (GeV)" 2 , R 2 = 3.18 (GeV)~ 2 , a' p = 0.25 (GeV)~ 2 , 
C= 1.5 , A = 0.139 . 
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Logi (E/GeV) 



FIG. 1: The three neutrino fluxes in the QGSM in the pp interactions as a function of the neutrino 
energy E. For the tau neutrino, the pQCD result is also presented for comparison. The three 
neutrino fluxes in the pp — * ti — * v + X, where t is the top quark, calculated in the pQCD, are also 
shown clS cl function of the neutrino energy E. 
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FIG. 2: The intrinsic astro-physical neutrino flavor ratios as a function of the neutrino energy E in 
the pp interactions. More details are given in the text. 
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FIG. 3: The ratios of the intrinsic neutrino fluxes, the Rp/ e and the R®/ e , for three different cosmic- 
ray flux spectrum exponents as a function of the neutrino energy E. The two intrinsic astrophysical 
neutrino flavor ratios are essentially independent of the cosmic-ray flux spectrum exponent. The 
8 = —2.75 case is also shown for comparison. More details are provided in the text. 
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FIG. 4: The oscillated astrophysical neutrino flavor ratios as a function of the neutrino energy E. 
For L(pc)/-E(GeV) > 10~ 10 , the three ratios enter into a region of the relatively rapid and small 
amplitude oscillations centered at the abscissa axis value 1. 
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FIG. 5: The distance to the astrophysical neutrino source L in units of parsecs as a function of 
the neutrino energy E in GeV, using the criterion L(pc)/£'(GeV) > 10 -10 . The shaded region 
indicates the L and the E value range for which the astrophysical neutrino flavor ratios may be 
1:1:1 originating in the pp interactions. 
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